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ABSTRACT 

Studies of quasar absorption lines reveal that the low density intergalactic medium at z ~ 3 is enriched 
to 10 -3 — 10 -2 solar metallicity. This enrichment may have occurred in an early generation of Population 
III stars at redshift z 10, by protogalaxies at 6 <; z <; 10, or by larger galaxies at 3 <; z < 6. This 
paper addresses the third possibility by calculating the enrichment of the IGM at z ^ 3 by galaxies 
of baryonic mass ;> 10 8 5 M Q . We use already completed cosmological simulations to which we add a 
prescription for chemical evolution and metal ejection by winds, assuming that the winds have properties 
similar to those observed in local starbursts and Lyman-break galaxies. Results are given for a number of 
representative models, and we also examine the properties of the galaxies responsible for the enrichment 
as well as the physical effects responsible for wind escape and propagation. We find that winds of 
velocity ;> 200 — 300 km s -1 are capable of enriching the IGM to the mean level observed, though many 
low-density regions would remain metal free. Calibrated by observations of Lyman-break galaxies, our 
calculations suggest that most galaxies at z > 3 should drive winds that can escape and propagate to 
large radii. The primary effect limiting the enrichment of low-density IG gas in our scenario is then the 
travel time from high- to low-density regions, implying that the metallicity of low-density gas is a strong 
function of redshift. 

Subject headings: cosmology: theory — intergalactic medium — galaxies: abundances, starburst 



1. INTRODUCTION 

The detailed comparison of quasar absorption spectra 
with the predictions of cosmological hydrodynamic sim- 
ulations has established that the Lya absorption 'for- 
est' is caused by a smoothly-fluctuating neutral compo- 
nent of the intergalactic medium (IGM) (e.g., Cen et al. 
1994; Zhang, Anninos, & Norman 1995; Hernquist et al. 
1996). The simulations also reveal a strong correlation 
between absorber H I column density and gas overdensity 

5 = p g as/(/Ogas), so that the study of AT (HI) < 10 15 cm~ 2 
absorbers gives information about the IGM for 5 <; 10 at 
z £ 3 (e.g. Zhang et al. 1998; Dave et al. 1999). Further- 
more, studies of CIV, Si IV and OVI have revealed that 
the low-density IGM is not - as was long expected - chemi- 
cally pristine, but has been enriched to ~ 10~ 25 Zq (with 
a scatter of - 1 dex) down to N(H I) » 10 14 5 (6 ~ 5) 
(e.g., Meyer & York 1987; Songaila & Cowie 1996), and 
perhaps to even lower overdensities (e.g., Cowie & Songaila 
1998; Ellison et al. 2000; Schaye et al. 2000). 

This enrichment may have occurred during a Popula- 
tion III star-formation epoch at z ^> 5 (e.g. Carr, Bond 

6 Arnett 1984; Ostriker & Gnedin 1996; Haiman & Loeb 
1997; Abel et al. 1998), later through the disruption of 
small 'protogalaxies' at 6 — 8 < z <; 10 — 20 (e.g. Gnedin 
& Ostriker 1997; Gnedin 1998; Madau, Ferrara & Rees 
2001), or even more recently by the ejection of metals from 
fairly massive galaxies at z <J 6 — 8. Earlier enrichment 
might allow more uniform IGM pollution (depending on 



how strongly biased the star formation is), but is con- 
strained by the amount of metal that can be generated 
at early times; conversely, sufficient metal is easily avail- 
able at late times, but the larger inter-galaxy spacing may 
preclude a homogeneous metal distribution (which obser- 
vations may or may not indicate). 

In this paper, we address enrichment of the IGM by 
supernova-driven winds from galaxies of (gas+star) masses 
,> 10 8 5 at z < 8. Our calculations employ a smoothed- 
particle hydrodynamics (SPH) simulation with 128 3 dark 
matter particles and 128 3 SPH particles in a (17Mpc) 3 
box, adopting a cosmology with fl\ = 0.6, = 0.047 and 
Q m = Qcdm + ^6 = 0.4. The simulation is described in 
more detail in Aguirre et al. (2001) and in Weinberg et al. 
(1999). These simulations do not include winds, but we use 
a parameterized prescription, described in § |, to add met- 
als in a way that models their ejection and distribution by 
supernova-driven galactic outflows. This post-processing 
method allows us to easily vary our assumptions and pa- 
rameters, and we give results for a number of representa- 
tive models. We then compare the resulting intergalactic 
metal distribution to the quasar absorption line observa- 
tions and discuss the implications for wind enrichment of 
the IGM. 

2. METHOD 

The method by which we calculate IGM enrichment is 
discussed in detail in Aguirre et al. (2001). Briefly, our 
method employs a limited number of outputs from already- 



2 



completed SPH simulations that include star formation. 
Each unit of forming stellar mass instantaneously gener- 
ates y* units of metal mass (y* = Zq is assumed for all 
models). We then deposit this metal mass in nearby gas 
particles as follows: 

1. Particles are grouped 1 into bound objects ('galax- 
ies'). The mass of star formation is tallied for 
a given galaxy, for which we also compute the 
mass, star-formation rate (SFR), and an area A — 
44(AT/10 n M©) - 76 kpc 2 based on an empirical rela- 
tion between disk area and mass (Gavazzi, Peirini & 
Boselli 1996). 

2. Galaxies with SFR/ A > SFR crit 

0.1 MQyr" 1 kpc" 2 are assumed to coherently drive 
steady-state, mass-conserving winds with velocity 
v out at the center of star formation and with a 
mass outflow rate of to = 1.3x(600kms" 1 /w ou t) 2 
times the SFR. Here \ parameterizes the frac- 
tion of supernova kinetic energy in the wind, if 
we assume 10 51 ergs are released per 100 M Q of star 
formation. Based on observations that generally 
to ~ 1 MQyr" 1 (Heckman et al. 2000), we use \ = 1. 
We allow galaxies with SFR/A < SFR crit to drive 
winds also, but with the wind energy attenuated by 
[(SFR/A)/ (SFR cr it)] 4 - 

3. We divide the volume about the center of star for- 
mation into N a regions of equal solid angle about 
directions (6^,^), with ~ 16 galactic gas particles 
per region (thus N a increases with galaxy mass). 

4. For each angle (8i,<fii) we integrate the equations of 
motion of a 'test shell' beginning at radius tq con- 
taining 10% of the galaxy's mass. The shell at radius 
r feels accelerations due to gravity, the ram pressure 
of the wind (attenuated by the potential difference 
between r and ro), the ambient thermal pressure, 
and the sweeping up of a fraction e cn t of the ambient 
medium (matter could be left behind by the shell 
either if the ambient medium is inhomogeneous, or 
as the shell fragments and leaves part of itself with 
lower velocity. We take a fiducial value of e cn t = 0.1). 
All accelerations are computed using the simulation 
particles contained in the solid angle about (9i,(f>i) 
(see Aguirre et al. 2001 for details). Eventually, the 
shell stalls (i.e. has a negligible velocity with respect 
to the local medium) at a radius r s t a ii(#, </>)■ We also 
stop the shell if/when its propagation time exceeds 
the time between its launch and the redshift at which 
we quote our results. 

5. In each angular region, metals are distributed among 
gas particles within r s t a ii(#, </>) such that the metal 
mass in a shell at r of thickness dr is cx (r/r sta ii) a dr. 
We generally take a — 3, but the resulting enrich- 
ment is rather insensitive to a for 2 <^ a <J 5 (Aguirre 
et al. 2001). 

6. If r sta ii < 2r for some angle, the corresponding 
metal mass is not distributed as per the wind pre- 
scription, but is deposited 'locally' over the 32 gas 

1 The grouping is done using the SKID package, publicly available 
2 The yield can be constrained using the metallicities of observed £ 



particles nearest to its progenitor star, using the SPH 
smoothing kernel (see, e.g., Hernquist & Katz 1989). 



The process is repeated for each galaxy, at each time 
step. New stars are formed with the metallicity of the gas 
from which they form. Metals accumulate in gas particles 
(which move between time steps), and the distribution of 
metals in stars and in galactic and intergalactic gas can be 
assessed at each time step. 

3. RESULTS 
3.1. Enrichment of the IGM 

We first give results for six illustrative models. The first 
four have all their parameters fixed (to the 'fiducial' val- 
ues given above) except for the wind outflow velocity iwt, 
which we take to be 600 ± 200kms _1 , 300 ± 200 km s" 1 , 
200±100kms _1 , or 100±50kms _1 (uniformly distributed 
in the given range for each model), independent of galaxy 
mass (as suggested by Heckman et al. 2000). Next, we fix 
^out = 300±200kms _1 and vary the entrainment fraction 
e ent from 1% to 100%. 

The key results of the calculations are shown in Fig. |[ 
The top two panels give a sparse sampling of individual 
particle metallicities, versus the gas overdensity 5. The 
stellar yield j/» is uncertain by perhaps a factor of two, 
and all of the curves could be scaled vertically for a higher 
assumed value. 2 The metallicity at <5 <; 10 4 could also be 
(roughly) scaled down by a factor (1 — Y lct ) if a fraction 
Yrct of the metal produced is assumed not to be incorpo- 
rated into the wind. Panel A shows that in the 600 km s" 1 
model, the enrichment is fairly uniform. A significant frac- 
tion of particles (shown in a bar at the bottom with an 
artificial 1 dex of scatter added) have zero metallicity, but 
the rest have a metallicity spread at a given overdensity 
of ^1-2 dex. In contrast, the 300 km s" 1 model (panel B) 
has some metal rich particles at low densities, but many 
more pristine regions that the winds have not reached. 

Panels C and D show the median and mean gas metal- 
licities vs. overdensity for all six models. The shaded box 
indicates the rough range of metallicity found in most ab- 
sorption systems for which metal lines are detected. The 
area with vertical lines indicates the region where the 
metallicity is agreed upon by independent groups (e.g., 
Rauch, Haehnelt, & Steinmetzet 1997; Dave et al. 1998). 
The horizontal lines indicate yet smaller overdensities to 
which this metallicity may extend (e.g., Ellison et al. 2000; 
Schaye et al. 2000). Even if the observational results were 
complete, neither the mean nor the median particle metal- 
licities we quote are directly comparable to the observa- 
tional results, but should bracket the quantity that is com- 
parable. The mean would be too high because (especially 
in cases with inhomogeneous distributions) it tends to be 
dominated by a small number of high- metallicity regions; 
the median gives information about the number of lines 
that would be metal rich, but is an underestimate because 
an absorption line would generally probe more than one 
simulation gas particle, decreasing the scatter. 

i jittp: / /www-hpec. astro. washington.edu/toole . 
laxies; see Aguirre et al. (2001). 
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Enrichment of the IGM at 
gas density. Panel A: Random subsample (1 in 500) of particle 
metallicities for wind model with v ou t = 600 km s -1 , versus over- 
density 5 = p g as/{Pgas). Top axis (here and in all panels) gives ap- 
proximate log N (HI), using the relation of Schaye (2001; assuming 
his fiducial parameters). The solid line shows the median metal- 
licity versus S. B: As for panel A, but in a v ou t = 300 km s -1 
model. C: Median metallicity versus <5 for wind models with 
Vout = 100, 200, 300, 600 km s -1 , and for 'minimal' and 'maximal' 
outflow models (as described in text). The shaded box roughly in- 
dicates the metallicity of low-column density Lya clouds. D: As 
for panel C, but mean metallicities are plotted. E: As for panel C, 
but gives mean metallicities times the fraction of baryons at a given 
5, showing the contribution by components with different <5 to the 
cosmic metal density. The thick solid line shows the distribution 
assuming constant metallicity (with the same total metal mass). 

In panel E the gas metallicity is multiplied by the frac- 
tion of gas at the given overdensity, showing what fraction 
of all cosmic metal lies at a given density. For the models 
that can plausibly account for the observed metal lines in 
the IGM, - 10 - 30% of all metal is extragalactic (EG), 
i.e. is not in cooled gas or stars in galaxies 

The curves show that the assumed wind velocity at in- 
ception is quite important in the resulting IGM enrich- 
ment. Slow winds tend to be confined, as can be seen 
from the small fraction (2%) of all metals outside of galax- 
ies in the 100 km s -1 model (see panel E). Those that 
escape tend to stall at relatively small distances where 

3 This assumes that the missing low-mass galaxies are distributed like the higher-mass galaxies. If low mass galaxies are strongly biased 
toward appearing in low-density regions, then the results may differ significantly when they are included. 



the density is still high, and fail to effectively pollute the 
low-density IGM. Increasing v out both increases the to- 
tal EG metal density (raising the EG fraction to 50% for 
Vout — 600 km s -1 ) and shifts the metals to lower den- 
sity regions (increasing the S = 1 mean enrichment by a 
factor of several hundred). Changes in the entrainment 
fraction also have a fairly large effect. If (instead of the 
fiducial 10%) only 1% of the ambient gas were entrained, 
20% of all metals would be outside of galaxies in our fidu- 
cial 300 km/s model; while only 4% would be extragalac- 
tic if all ambient gas were entrained. Changing our other 
model parameters within a range of reasonable values has 
a smaller but still substantial effect; doubling SFR cr ;t or 
changing \ by a factor of two changes the ejection frac- 
tion by <J 30% and changes the mean enrichment at 5 = 1 
by <, 50% (Lowering SFR cr ;t would have little effect since 
most galaxies at z > 3 are already driving winds for our 
assumed value.) 

Two methodological points concerning our results merit 
discussion. First, resolution tests reveal that galaxies with 
;> 6 x 10 8 M Q in baryons are well-resolved. Galaxies be- 
low this limit exist (and are used in the analysis) but will 
be under-represented both in their number and in their 
star-formation rates (Weinberg et al. 1999). Studies of 
dwarf galaxies at high z have concluded that at z <^ 6, 
galaxies of baryon mass ^5x 10 7 M are suppressed by 
photo-evaporation (e.g., Ferrara & Tolstoy 1999), leaving 
a (redshift-dependent) 'gap' in mass where galaxies may 
exist but are not resolved by our simulations. 

Second, our method necessarily neglects the effect of the 
outflows on the galaxies. Galaxy formation studies with 
ineffective feedback generally predict more low-mass galax- 
ies than are observed, and the simulations used here may 
over-estimate the cosmic SFR and stellar mass for this rea- 
son (see Weinberg et al. 1999 and Aguirre et al. 2001). 
While we cannot quantify exactly the effect these two (op- 
posing) inaccuracies will have on our results, comparison 
of the mass functions for different numerical resolutions 
indicates that we can probably hope to estimate the en- 
richment by z <; 6 galaxies to within a factor of a few. 3 

Because IGM enrichment is so sensitive to v ou t, we can 
set a lower limit on that quantity if the observed enrich- 
ment is primarily due to winds of the type discussed in 
this paper. To achieve the mean metallicity of 10~ 2 ' 5 Zq 
observed in JV(HI) ~ IO 14 5 absorbers, w ou t ^ 200 km s" 1 
is necessary, and w ou t ^ 300 km s -1 would be needed to 
enrich absorbers down to iV(HL) ~ 10 135 with a similar 
metallicity (though this may over-enrich higher column 
density systems). A higher velocity is probably necessary 
if the observed enrichment is fairly uniform; for a median 



metallicity of - 10 5 Z Q down to AT(HI) ~ 10 5 cm" 2 , 
a wind speed of v out ;> 300 km s -1 is probably required, 
though a more detailed analysis is necessary for a robust 
comparison. These conclusions depend on the assumed en- 
trainment fraction (and more weakly on the other parame- 
ters), but changing these assumptions cannot increase the 
enrichment very much because enrichment is always lim- 
ited by the available time between shell launch and z = 3. 
Thus if the observed enrichment turns out to be uniform at 
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column densities AT(HI) <C 10 5 cm 2 , enrichment from 
smaller and/or higher redshift galaxies would be required. 



3.2. Physics of outflows 

In addition to providing estimates of intergalactic en- 
richment, our calculations also shed light on the physics 
of wind escape. Figure ^ gives information about the out- 
flows propagating between z ~ 4 and z = 3 from indi- 
vidual galaxies in the v out — 300 km s" 1 model. Panel 
A shows the shell stalling radius for each shell about a 
sample of 500 simulation galaxies; each point corresponds 
to one propagation direction and each vertical 'stripe' to 
a galaxy. At 3 < z < 4 about half of all galaxies are 
driving winds that reach r sta ii lOkpc. The variation 
in r s t a ii for a given galaxy is due to anisotropics in the 
galaxy and IGM; variations in r sta n among galaxies of a 
given mass are due to combined variations in the outflow 
velocity, the SFR and the properties of the nearby IGM. 
Since we have assumed that wind velocity is independent 
of galaxy mass, the most massive galaxies drive winds in- 
effectively because of their relatively large potential wells. 
The cutoff at r s taii ~ 100 kpc is due to the finite time 
between z ~ 4 and z = 3 (where we stop each shell's in- 
tegration). This is shown in Panel B, which displays the 
stalling time for the shell segments (one point for each an- 
gle). As shown in Panel C, at z — 4 almost all of the 
galaxies exceed the assumed critical SFR/(area). But as 
long as this threshold is satisfied, the areal SFR does not 
strongly influence how well winds escape. 

Panel D allows a comparison of our results for wind 
escape with two other investigations. The points in the 
figure indicate the supernova energy injection rate for 
each galaxy in units of 10 38 ergs, converted from star for- 
mation rates by assuming that one supernova releasing 
10 51 i?5i ergs forms per per lOOMioo M of star formation. 
The vertically shaded region indicates the range of SFRs 
measured in bright, wind-driving Lyman-break galaxies by 
Pettini et al. (2001), converted in the same way. The 
solid line is a fit that extrapolates the simulation results 
to lower-mass galaxies, and the dotted line represents our 
fiducial critical SFR/area (i.e. we assume that galax- 
ies above the dotted line are driving winds coherently - 
though the winds may ultimately be retained as shown in 
panel A). Mac Low & Ferrara (1999; see also Ferrara & 
Tolstoy 2000) have studied galactic winds in dwarf galax- 
ies using numerical simulations, for the range of galaxy 
masses and supernova energy injection rates shown by the 
shaded box. They find that only the lowest-mass galaxies 
(10 6 M Q ) can 'blow-away' the bulk of their ISM, though 
higher mass galaxies can allow 'blowout' where hot, metal 
rich gas can escape. Their results agree roughly with the 
estimates of Silich & Tenorio-Tagle (2001), who find that 
winds can only escape spherical galaxies (corresponding to 
'blowout') for energy injection rates above the dashed line 
of Panel D; winds can escape flattened galaxies for rates 
above the dot-dashed line. According to our calculations, 
winds can escape the simulation galaxies (which are not 
highly flattened), even for entrainment fraction unity, from 
galaxies with mass 10 6 M . This would seem to conflict 
with the findings of Mac Low & Ferrara (1999), but the 
disagreement is probably due to the rather small energy 
injection rates they assume for their largest galaxies - 10 4 
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Fig. 3. — Enrichment of the IGM at z = 3 for the fiducial model 
with v = 300 km s -1 , with various effects turned off. A: Mean 
metallicities versus S for wind models with gravity, wind ram pres- 
sure, IGM thermal pressure, IGM infall, or the time limit turned 
off. B: As in panel A, but median metallicities. C-F: Wind stalling 
radius for wind propagating between z = 4 and z = 3, versus galaxy 
baryon mass, with various effects turned off (see text). 

times smaller than the rate indicated by the cosmologi- 
cal simulations. The simulated galaxy SFRs at high z are 
high enough to fulfill the (extrapolated) 'blow-away' crite- 
rion of Silich & Tenorio-Tagle for all galaxy masses, as are 
the SFRs of the Lyman-break galaxies (unless the galaxies 
have masses > 10 11 M .) 

In Figure |we show plots of both IGM enrichment and 
stalling radius vs. mass in the fiducial model with each 
effect on shell propagation turned off in turn. First (panel 
D) we remove the time constraint. This shows that some 
shell fragments from z = 4 galaxies would ultimately 
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the strip just below 10 9 yr corresponds to the time between z = 4 and z = 3. C: r sta n vs. total SFR/area of galaxy. D: Supernova energy 
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of star formation. Points show the simulated galaxies, to which the solid line is fit. The dotted line corresponds to our fiducial value of 
SFR cr ; t . The upper (vertically) shaded region is the range of SFRs (10 — 50MQ/yr) found by Pettini et al. (2001) converted in the same 
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propagate to ~ 1 Mpc (the cutoff comes from the distance 
d at which Hd = v out ). Panel E shows the results if the 
wind is turned off 10 8 years after the shell launch. This 
is roughly the minimal interval Td over which a starburst 
'event' could drive the shell, since the energy deposition 
rate is roughly constant for r e ~4x 10 7 yr after a coeval 
starburst (Tenorio-Tagle et al. 1999), and the wind (or 
thermal pressure) will still drive the shell for the time r c it 
takes the wind (or sound) to cross the shell radius. If the 
shell speed is comparable to the wind speed, r e ~ r c so 
Td ~ 2r e . The effect of truncating the wind is noticeable 
but slight. In panel F we have turned off gravity, and es- 
sentially all of the shells escape, showing that in this case 
gravity is the dominant effect in confining winds. Of the 
shells that still do not escape, about half are retained due 
to thermal pressure, and about half due to the sweeping up 
of matter. Panel G shows the effect of very low mass en- 
trainment: winds escape to slightly larger radii on average, 
but can actually increase the escape radius in some cases 
because higher entrainment adds more momentum to the 
shell. Entrainment actually causes as much deceleration 
on the shell as gravity; but since gravity also attenuates 
the wind (decreasing the acceleration of the shell), it has 
a larger overall effect. Panel H shows that taking into 
account infall from the IGM is not very important, but 
tends to slightly decrease the stalling radius. In summary, 
we find that only the hot wind, the entrainment of am- 
bient material and gravity are important in determining 
shell propagation in most cases at z = 3. More hot gas 
(especially in clusters) and stronger infall, however, might 
make thermal pressure and IG infall more important at 
low z. 

The most important factor limiting the enrichment of 
very low density regions is the time constraint as shown 
in Fig. [|. If we assume that the winds can propagate 'in- 
stantaneously' into the IGM, they reach quite low-density 
regions, and the enrichment at z = 3 is only slightly higher 
and more uniform than at z — 4. But if the time limit is 
respected, the enrichment spreads to significantly lower- 
density regions between z = 4 and z — 3 due to the extra 
Gyr of propagation time. 

4. DISCUSSION AND IMPLICATIONS 

Metal enrichment of the z ~ 3 IGM has also been in- 
vestigated by Gnedin & Ostriker (1997), Gnedin (1998) 
and Cen & Ostriker (1999) using simulations that include 



chemical evolution. Gnedin finds that supernova blow- 
apart of protogalaxies is ineffective, but that dynamical 
removal of metals can explain the metals in low-density 
regions - though this requires a metal formation rate at 
z ;> 4 that is about 10-50 times that observationally in- 
ferred at z ~ 3 — 4 (e.g., Stcidel et al. 1999) and implies 
more efficient dynamical removal of metals than found by 
Aguirre et al. (2001). 

Cen & Ostriker (1999) assume a fixed metal ejection 
fraction into regions of ~ 400 kpc (the resolution of their 
grid-based hydrodynamical treatment; see Cen & Ostriker 
2000). Their results for mean metallicity are similar to 
our v = 200 km s -1 model, though their metallicity distri- 
bution at fixed 5 shows almost no scatter at z = 3 (but 
becomes more inhomogeneous with time). In contrast, our 
calculations find a highly inhomogeneous metallicity dis- 
tribution, which becomes somewhat more uniform with 
time. This may be due to the rather low spatial resolution 
of their simulations or may indicate that stars are form- 
ing (and expelling metals into the IGM) in very low den- 
sity regions in Cen & Ostriker's simulations, in contrast to 
ours. Comparison of quasar absorption spectra to simu- 
lated spectra shows that ^0.5 dex of scatter in gas metal- 
licity (among absorbers in which metals are detected) may 
be necessary to explain the variation in observed CIV/HI 
ratios (Rauch et al. 1997; Hellstcn et al. 1997; Dave et al. 
1998), so this point merits further investigation, but will 
probably require the generation of synthetic spectra from 
our simulations. 

An important physical effect that is not addressed by 
our treatment is that of the outflowing winds on the IGM: 
if metals are transported directly from galaxies to very low 
density regions by movement of large quantities of gas, 
then this kinetic energy must be absorbed by the IGM. 
Yet the widths of low-column density Lya forest lines are 
typically 20 — 30 km s _1 , and line ratios indicate photoion- 
ization, not collisional ionization (e.g. Songaila & Cowie 
1996). Our method does not allow us to assess the effects 
of winds on the IGM (see Theuns, Mo & Schaye [2000] 
for some discussion), but we can make some general com- 
ments. If expanding shells of swept-up gas maintained 
their integrity and propagated to roughly the inter-galaxy 
spacing, they would completely re-make the structure of 
the Lya forest, except perhaps in deep voids. But this is 
not the expected scenario: upon leaving a galaxy, the wind 
and 'shell' will likely mix into a complex, multi-phase out- 
flow consisting of fragmented shell(s), the hot wind, dense 
clouds of entrained (but not swept-up) ambient medium, 
etc. Small, dense blobs might maintain fairly low temper- 
ature through efficient cooling, and could reach large radii 
before stalling and mixing with the IGM (and by defini- 
tion when the wind stalls it has small local velocity and 
will not greatly heat the medium). In addition, the hot 
wind itself will not be observable in absorption, and will 
flow more more easily in low-density directions. Hence it 
may not seriously disrupt the relatively dense filaments 
that give rise to much of the Lya absorption. 

These issues merit further study because our results, 
combined with observations of high-z galaxies, indicate 
that significant outflow from z > 3 galaxies seems in- 
evitable. Essentially all of the simulated galaxies at z > 3 
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Fig. 5. — Maximum stalling radius in each galaxy vs. star for- 
mation rate, at z = 4, for v ou t = 750 ± 500 km s -1 and entrainment 
fraction unity. The Lyman-break galaxies of Pettini et al. (2001) 
show SFRs of 10 - 50 M yr" 1 . 

display specific time-averaged SFRs exceeding the ob- 
served threshold for driving superwinds in local (starburst) 
galaxies - and stochastic star formation should result in 
much higher rates for short periods. Moreover, available 
observations of z ~ 3 galaxies reveal that they are driv- 
ing winds, with characteristic velocities exceeding those of 
our fiducial model: Pettini et al. (2001) find that Lyman 
break galaxies at z ~ 3 with SFRs of ~ 10—50 M /yr drive 
winds at ~ 100 — 800 km s -1 (as measured by absorption 
lines) or ~ 250 — 1250 km s -1 (in Lya emission). These 
galaxies correspond to the highest (~ 10 10 M Q ) mass sim- 
ulated galaxies (see Fig. || and Dave et al. 1999). The 
outflow signatures in the spectra are quite similar to those 
in the nearby wind-driving starbursts upon which we base 
our method, so the observed velocities should roughly cor- 
respond to our assumed initial outflow velocities. In this 
case, as shown in Fig. ^, the outflows tend to escape the 
galaxies even for an entrainment fraction unity, and prop- 
agate to large distances - none of the forces we include 
in our calculations would be capable of containing the 
outflow. Thus we predict that winds from Lyman-break 
galaxies do escape those galaxies' potentials and enrich 
the IGM, 4 and expect that this would hold for lower mass 
galaxies as well. 

5. CONCLUSIONS 

We have simulated the enrichment of the z = 3 IGM by 
galactic winds from ;> 10 85 M Q galaxies, assuming that 
most drive winds with properties similar to those observed 
in local starbursts and in Lyman-break galaxies. We find 
that for winds with outflow velocity v out > lOOkms -1 en- 
richment is significant, and that if w out ;> 300 kms , the 
winds can (roughly) account for the metallicity of the low- 
column density (AT (HI) ~ 10 14 5 cm -2 ) Lya forest. Higher 
velocities could enrich yet lower-density components of the 
IGM. 

Our treatment of wind-driven outflows includes forces 
due to gravity, thermal pressure, ram pressure of the am- 

4 If the Lyman break galaxies actually represent lower-mass, very strong 
it may then be less applicable to the rest of the z ~ 3 population. 



bient medium, and the hot wind itself. We find that at 
z ~ 3 thermal pressure and the the infall velocity of the 
IGM are generally negligible, but that outflows can be 
confined by gravity or ambient material for sufficiently 
high mass galaxies and/or low velocity winds. Our results 
suggest that the winds observed in Lyman break galaxies 
should escape and travel to fairly large distances. If the 
Lyman break galaxies represent high-mass, non-starburst 
galaxies, then the lower-mass z = 3 galaxies should drive 
outflows to large distances as well, unless outflow veloc- 
ity depends strongly on galaxy mass. The effect of such 
winds on the physical and thermal structure of the IGM is 
essentially unknown. The good agreement between simu- 
lations like these and the observed statistical properties of 
the Lya forest suggest that the effect cannot be extremely 
large. 

We find that the main effect preventing the enrichment 
of very low-density regions of the IGM is the time available 
for metal-rich material to propagate from galaxies (which 
are concentrated in filaments) into the voids. This leads 
to a picture in which a given set of galaxies enriches a pro- 
gressively larger cosmic volume with decreasing redshift, 
and we predict that if low-density regions are enriched by 
winds from ;> 10 85 M Q galaxies at z <J 6, the median 
metallicity of low-density regions should be a very strong 
function of redshift. 

The scenario described here, enrichment of the IGM by 
massive (;> 10 8 5 Mq) galaxies at modest (z <^ 6) redshift, 
is but one of several interpretations for the origin of metal 
in the Lya forest. Other possibilities, distinguished by 
epoch and mass scale, include enrichment by population 
III stars at z ~ 20, or the disruption of low-mass proto- 
galaxies at z ~ 10. How are we to decide which, if any, of 
these three scenarios is most relevant to the real Universe? 

Based on our modeling, we anticipate that these mech- 
anisms will predict differences in the metal distribution of 
the IGM that can be tested observationally. In addition 
to the mean enrichment, the dispersion in metallicity at a 
given overdensity and the trend of metallicity with over- 
density can both be measured. It seems natural to expect 
that metal deposition at relatively higher redshift will lead 
to a rather different distribution of metals than that in- 
ferred from our analysis. Thus, future modeling, combined 
with observational samples of increasingly large size, may 
discriminate between the various scenarios. To the extent 
that this goal can be realized, the metals discovered in the 
z r~j 3 Lya forest may provide a direct probe of the condi- 
tions responsible for the production of the first metals in 
the Universe. 
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ly star-bursting galaxies this conclusion is only strengthened, though 
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